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ABSTRACT 


Experimental measurements are made with a novel two wire thermocouple. 
Signals from two wires of unequal diameters are recorded from the thermocouple 
suspended in constant flow with a periodic temperature fluctuation. It is 
demonstrated that the reconstructed signal from the two wire thermocouple 
requires no compensation for to <. 2 0)1 where coi is the natural frequency of the 
smaller wire. A compensation factor is recommended for larger frequencies co > 2o>i. 

Theory and experimental measurements are compared with a novel three 
wire thermocouple. Signals from three wires of unequal diameters are recorded 
from the thermocouple suspended in constant flow with a periodic temperature 
fluctuation. It is demonstrated that the reconstructed signal from the three wire 
thermocouple requires no compensation for co<5coi where (Oi is the natural 
frequency of the smaller wire. The latter result represents a significant 
improvement compared to previous work with two wire thermocouples. A 
correction factor has also been derived to account for wires of arbitrary diameter. 

Measurements are recorded for multiwire thermocouples consisting of either 
two or three wires of unequal diameters. Signals from the multiwire probe are 
recorded for a reversing gas flow with a periodic temperature fluctuation. It is 
demonstrated that the reconstructed signal from the multiwire thermocouple 
requires no compensation provided co/coi < 2.3 for two wires or (o/coi < 3.6 for three 
wires where (Oi is the natural frequency of the smaller wire based on the maximum 
gas velocity. The latter results were possible provided Fourier transformed data 
were used and knowledge of the gas velocity is available. 
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TWO WIRE THERMOCOUPLE: FREQUENCY 
RESPONSE IN CONSTANT FLOW 



1. INTRODUCTION 


The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine. Historically, these temperatures have 
been measured with thermocouples. 1 The advantages of thermocouples are their 
low cost, reliability, and simplicity since they do not require optical access or 
elaborate support electronics. However, the design of a thermocouple represents a 
compromise between accuracy, ruggedness, and rapidity of response. 

For example, the measurement of fluctuating temperatures in the high-speed 
exhaust of a gas turbine engine combustor is required to characterize the local gas 
density gradients or convective heat transfer. 2 Although thermocouples are suitable 
for the measurement of high-frequency temperature fluctuations (<1 kHz) in a 
flowing gas or liquid, the measured signal must be compensated since the frequency 
of the time-dependent fluid temperature is normally much higher than the natural 
frequency of the thermocouple probe. 3 Moreover, use of a single wire thermocouple 
in constant velocity flows requires knowledge of the fluid velocity and properties 
(e.g., viscosity, density, etc.) to determine the natural frequency. 

The present paper describes the performance of a novel two wire 
thermocouple of unequal diameters that does not require compensation at lower 
fluid temperature frequencies nor any knowledge of the fluid velocity or properties. 
The results of experimental measurements are presented along with the suggested 
procedure for the reduction of the data from the two wire thermocouple as shown 
in fig. 1. 
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2. ONE WIRE 


Use of a single wire thermocouple requires knowledge of the fluid velocity 
and properties to determine the natural frequency. The latter quantity is necessary 
to establish a frequency dependent compensation factor for the measured signal. 


2.1 Steady State Response 

The conservation of energy for a single wire is : 


% = “.A- T > 


( 1 ) 


where the natural frequency 


Here: 


g 
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T 
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4h 
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(2) 


= natural frequency (sec -1 ) 

= temperature of the thermocouple ( °K) 

= temperature of the gas (°K) 

= time (sec) 

= wire density (k g /m 3 ) 

= heat capacity of wire (J/K g - "K) 

= heat transfer coefficient (W/m 2 - ‘K) 

= wire diameter (m) 

It is convenient to rewrite the heat transfer coefficient in terms of the Nusselt 
number 


h = 


— Nu 
D 


(3) 
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where 


Nu = 


thermal conductivity of the gas (W/m - °K) 
Nusselt number. 


Thus, the natural frequency in Eq. (2) becomes 


4k Nu 

(0 = — ^ , 

(pc)D 2 


(4) 


where 


Nu = C Re m Pr 


m 0*^/3 


(5) 


and 


Re = 


YD 

\)„ 


Pr = 


^ c p 


We now assume a periodic variation in the temperature of the form 


T = p itot 

g 


(6) 


and a wire temperature 

T = Ae i0 * • (7) 

Substituting Eqs. (6) and (7) into Eq. (1), one obtains a wire temperature or first order 
response to a simple harmonic gas temperature fluctuation of the form 
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i(cot+<|> ) 
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where the phase angle 

0 n = - tan -1 (to/ co n ). 


(9) 


2.2 Error 

The natural frequency of a single wire from Eqs. (4) and (5) can be written in 
the form. 


m r>~l / 3 


CQ„ = 


k g C 0 4Re m Pr 


(pc)D 2 


( 10 ) 


or separating the wire diameter 


m-2 


co„ = CD 


where the coefficient C depends on both fluid and wire parameters or 


(ID 


4k C V m Pr 1/3 
C = - 8 ° 


(pc)D 


( 12 ) 


Substituting values corresponding to a type K thermocouple and a gas velocity V = 
10 m/s for air at standard conditions into Eqs. (11) and (12), one obtains values of the 
natural frequency: 

(O n (sec 1 ) D (m) 

33.8 50x1 O' 6 

75xl0- 6 


18.3 
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8.5 125xl0- 3 * * 6 

The values of the natural frequency to n depends on the choice of the 


exponent m in the empirical expressions Eqs. (4) and (5). Since a variation can exist 
in the value of m 4 , a computation has been made to determine the effect of the 
magnitude of m on the amplitude of the steady state response. Thus, from Eq. (8) 
one computes an amplitude ratio defined as 


T(m) 

T(.5) 


-IV 2 


(0 


co n (.5) 


J 


f 


(0 


+ 1 


+ 1 


(13) 


where co n (m) is determined from Eqs. (11) and (12) for a wire diameter D = 50 x 10' 6 
m (50pm). 

Figure 2 is a plot of Eq. (13) where it is demonstrated that a significant error 
can occur in the computation of the amplitude. In particular, an error of 40% in the 
amplitude will occur with a change of 20% in the magnitude of the exponent m. 
The problem can be significant for angular frequencies to > a) n 


3. TWO WIRE 


The use of a single wire thermocouple requires knowledge of the fluid 

velocity and additional properties such as viscosity or density. There are three 
unknowns that appear in Eqs. (1) and (11): a) the coefficient C given by Eq. (12) b) m 

the exponent of the wire Reynolds number in Eqs. (10), (11) and (12) and c) the gas 
temperature T g . 
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3.1 Determination of Gas Temperature 

Consider two thermocouples of unequal diameters Di <D 2 as shown in fig. 1. 
The conservation of energy equation for both reduce to 

= CD™ -2 (T.-T,) 04) 

dt 8 

^ = CDr 2 (T.-T,)- 05) 

dt 8 

Eliminating C and solving for T g/ one obtains 



where Ti, T 2 are the measured temperatures and their derivatives from wires of 
diameter Di, D 2 , respectively. Use of Eq. (16) only requires an estimate of the 
exponent m and knowledge of the wire diameters Di and D 2 . 

The construction of the gas temperature from Eq. (16) also requires knowledge 
of the derivatives of the wire temperatures. Calculations of these derivatives may 
amplify noise in the sampled data. Since Eqs. (14) and (15) are linear equations, we 
choose to use transformed data profiles. Taking the Fast Fourier Transform (FFT) 
of Eqs. (14) and (15), one obtains 


icoT, = CD”- ! (T - T,) 

(17) 

i<oT7 = CD"- j (T s - T 2 ) 

(18) 


or solving for T g 
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T k = T, 


1 - 


V D 2 ) 


m-2 




(Vj) 

— 



<^ 2 > 


in -2 


(19) 


The inverse transform of T g is the reconstructed gas temperature or 


t b = fft-'[t 8 ], (20) 

3.2 Error 

Assuming the wire diameters are known accurately, the choice of the 
exponent m will affect the accuracy of the prediction of the gas temperature T g as 
shown in Eq. (16). The effect of the exponent m on the predicted amplitude of T g is 
shown in fig. 3 which plots the amplitude ratio T g (m)/T g (.5) as a function of the 
normalized angular frequency to/toi . In fig. 3 values of Di = 50pm and D 2 = 75pm 
were chosen including a value of the exponent m = 0.5. Substitution into Eq. (8) 
provided values of Ti(Di = 50p) and T 2 (D 2 = 75pm) in Eq. (16). 

As indicated in fig. 3, improper choice of the exponent m will create 
significant error in the prediction of the gas temperature. The error is magnified at 
larger frequencies to > 3 coi for ~ 25% error in m while the onset of the error occurs at 
to > 10 coi for a smaller ~ 10% error in m. Compared with similar computations for a 
single wire thermocouple in fig. 2: a) the error is much smaller b) the error for the 
two wire probe is shifted to larger frequencies to > toi and c) the error is reasonably 
constant for all choices of the exponent m. 

Similar computations were made for a range of wire diameter ratios 1/2 < 
D 1 /D 2 <_9/10 with Di = 50pm. As shown in fig. 4, the error decreased by ~ 12% with 
wire diameters nearly equal. 
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4. EXPERIMENT 


4.1 Apparatus 

In the present experiment, thermocouple sensors are exposed to a 
constant velocity air stream (<19 m/s) of varying temperature. In particular, the 
dynamic response of the thermocouple is measured for a periodic temperature 
profile of varying frequency. A rotating wheel configuration is used to deliver the 
test air stream to the proposed sensors. A similar experimental apparatus was 
described in detail by Elmore et al 5 and Forney et al. 6 

A schematic of the rotating wheel apparatus used in the present experiment is 
shown in Fig. 5. As the wheel rotates, holes pass the two air supply tubes (3/4-in. i.d. 
copper) that allow slugs of hot (~55°C) and cold (~30°C) air to alternately enter a 
transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing a periodic temperature profile covering a range of gas temperature 
frequencies from roughly 1 to 30 Hz. A 3/4-1 /2 in. smooth copper adapter was 
inserted at the end of the coalescing tube (location of thermocouple) in Fig. 5 to 
increase the air velocity to 18 m/s. 

The analog temperature signal is digitized with a Data Translation DT-2801 
A/D board mounted in an expansion slot of an IBM AT compatible computer as 
shown in Fig. 5. ASYST software was loaded onto the hard disk of the personal 
computer and this provided a flexible system for data storage, manipulation, and 
display. The true temperature profile of the airstream is measured with a constant 
current anemometer (TSI 1054- A) and sensor (1210-T1.5). 

The two type K thermocouple wires were of diameter Di = 50pm and D 2 = 
75pm as shown in fig. 1. The thermocouple junctions were fabricated by the 
Research Instrumentation Branch at the NASA Lewis Center. The wires were cut 
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with a razor blade to produce a flat edge perpendicular to the wire axis. The wire 
segments were then mounted in a fixture such that the two faces of the junction 
were held together by springs. The junction was then formed by laser heating. The 
laser used was a KORAD model KWD Nd: YAG laser, operating at a wave length of 
1060 nm. The power settings depend on the wire size and material, but for the 
thermocouples described in this report, a pulse duration of 4 ms was used delivering 
a total energy of approximately 2J. 

4.2 Procedure 

To calibrate the two thermocouples and anemometer, the air velocity at the 
exit of the coalescing tube in Fig. 5 (location of thermocouple) was maintained at 
roughly 17.5-18 m/s. The voltage regulator was adjusted such that the temperature 
of the hot line was maintained at an elevated temperature of ~55°C. One hole of 
the rotating wheel was positioned in a stationary position over the unheated line. 
The voltage output from both the thermocouple and constant current anemometer 
were recorded at the input terminals to the A/D board. The true temperature was 
also recorded with a thermometer suspended in the high-velocity air stream. This 
procedure was repeated for several temperatures with the hole positioned over the 
hot line covering an air temperature range from 30 to 55 'C. It was found that both 
type K thermocouples had a constant calibration factor CC/V) over the small 
temperature range. The constant current anemometer, however, provided a 
nonlinear response with the calibration factor varying with air temperature. Thus, 
a least squares curve fit was determined to provide the gas temperature with the 
anemometer output. The calibration factors described above were then installed on 
the ASYST data acquisition software. 

The thermocouple and constant current anemometer are mounted in the 
constant velocity air stream with both sensing elements parallel and separated by 


12 



approximately 1 mm. The period or frequency of the periodic temperature profile is 
controlled over the range from 1.5 to 30 Hz with adjustment of the motor connected 
to the rotating wheel shaft. 

ASYST software was developed to acquire temperature data sequentially from 
the thermocouples and constant current anemometer. The data are digitized for 
three channels at a sampling rate of 1024 Hz per channel for a total sample time of 1 
s. After data acquisition, the Fast Fourier Transform (FFT) is taken for each channel. 
The largest peak (first harmonic) in the amplitude of the FFT is located for each 
channel which provides the frequency of the temperature profile. The ASYST 
software also records the amplitude ratio of the first harmonic for the 50|im 
thermocouple - to - anemometer output. The amplitude ratio was recorded for each 
frequency setting of the rotating wheel and provided a value for the natural 
frequency coi of the small wire. 

The FFT from each thermocouple channel was substituted into Eq. (19) along 
with values for the diameter ratios and an estimate of the exponent m by the ASYST 
code. The inverse transform of T g by the ASYST code provided a reconstruction of 
the true gas temperature. 

5. RESULTS AND DISCUSSION 

Temperature profiles from the two thermocouples and anemometer were 
recorded for several angular frequencies of the rotating wheel. Figure 6 illustrates 
three profiles covering a time span of 0.3 s for a wheel angular frequency of co = 74 s' 1 
and a natural frequency coi = 42 s* 1 for the smaller 50|im wire. Here, the attenuation 
and phase shift are evident for the two interior profiles corresponding to the two 
wire thermocouple response. 
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The FFT of the two thermocouple signals was determined and the results 
were combined as suggested by Eq. (19). The inverse FFT was taken of T g and the 
result was plotted along with the anemometer output as shown in fig. 7. As 
indicated, the anemometer output and the reconstructed signal from the two 
thermocouples are nearly equal. 

A series of experiments were performed with the wheel frequency covering a 
range of values from 2.5 to 25 Hz. The amplitude ratio of the reconstructed signal to 
the gas temperature (anemometer output) was plotted versus co/a>i as shown in fig. 
8. Each data point in fig. 8 represents the average ratio from two consecutive peaks 
derived from plots similar to fig. 7. 

For fixed exponent m = 0.4, fig. 8 illustrates that the reconstructed signal 
requires no compensation until to - 2(H\. For larger frequencies to > 2coi, the data 
conform to the empirical function for the amplitude ratio 


FFT ' 1 



|T, 



1.36 - 0 . 265 ( 0 )/©,) 


( 21 ) 


where the numerator is derived from Eqs. (19) and (20) and the denominator is the 
output amplitude from the anemometer. Also shown for comparison in fig 8 is the 
first order response for a single thermocouple wire. 

It is possible to predict the natural frequency of the small wire coi from Eqs. 11, 
17 and 18. Eliminating the transform of the gas temperature Tg from Eqs. 17 and 
18, one obtains 


( 0 , = ICO 


(D, /P 2 ) m ~ 2 T2-T1 
(T1-T2) 


( 22 ) 
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Equation (22) is a useful expression provided the first harmonic frequency from the 
transform of the gas temperature Tg is below 2toi. Evidence to support the latter 
must be derived from prior knowledge and estimates of the gas velocities, properties 
and the temperature fluctuation frequency. Knowledge of the natural frequency coi 
from Eq. (22) would allow one to use the compensation factor Eq. (21) for the 
amplitudes of the higher harmonics of Tg where co > 2coj. 

From Eq. (19), the amplitude ratio should depend on the exponent m. 
Additional data were recorded for m = 0.3, 0.4, 0.5 and 0.6 and the results are shown 
in fig. 9. It is apparent that m has little influence on the amplitude ratio. 
Theoretical results suggest that there must be a single value of m that would 
provide an amplitude ratio of one for all frequencies. However, the empirical 
nature of Eq. (5) for the Nusselt number and sampling error provide the attenuation 
in Fig. 8. 

6. CONCLUSIONS 

Experimental measurements are made with a novel two wire thermocouple. 
Signals from two wires of unequal diameters are recorded from the thermocouple 
suspended in constant flow with a periodic temperature fluctuation. It is 
demonstrated that the reconstructed signal from the two wire thermocouple 
requires no compensation for to <. 2 o)i where toi is the natural frequency of the 
smaller wire. A compensation factor is recommended for larger frequencies co > 2coi. 
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8. FIGURES 





Fig. 1 - Schematic of two wire thermocouple. 
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Fig. 2 - Effect of exponent m on the amplitude for single wire. 
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Fig. 3 - Effect of exponent m on the amplitude for a two wire 
thermocouple. 
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Fig. 4 - Effect of wire diameter ratio on the amplitude. 
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Fig. 5 - Schematic of rotating wheel apparatus. 


22 




xE0 



23 


Time (sec) 

Fig. 6 - Temperature profiles for the gas (large amplitude) and two type 
K thermocouple wires. Two wires: Di = 50pm, (Oi = 42 s* 1 , D 2 = 
75pm. Gas temperature frequency: to = 74 s _1 . 
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Fig. 7 - 
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Fig. 8 - Amplitude ratio of two wire thermocouple output-to-true gas 
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Fig. 9 - Effect of exponen 
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1. INTRODUCTION 


The evaluation of jet engine performance and fundamental studies of 
combustion phenomena depend on the measurement of turbulent fluctuating 
temperatures of the gas within the engine. Historically, these temperatures have 
been measured with thermocouples. 1 The advantages of thermocouples are their 
low cost, reliability, and simplicity since they do not require optical access or 
elaborate support electronics. However, the design of a thermocouple represents a 
compromise between accuracy, ruggedness, and rapidity of response. 

For example, the measurement of fluctuating temperatures in the high-speed 
exhaust of a gas turbine engine combustor is required to characterize the local gas 
density gradients or convective heat transfer. 2 Although thermocouples are suitable 
for the measurement of high-frequency temperature fluctuations (<1 kHz) in a 
flowing gas or liquid, the measured signal must be compensated since the frequency 
of the time-dependent fluid temperature is normally much higher than the natural 
frequency of the thermocouple probe. 3 Moreover, use of a single wire thermocouple 
in constant velocity flows requires knowledge of the fluid velocity and properties 
(e.g., viscosity, density, etc.) to determine the natural frequency. 

The present paper describes the performance of a novel three wire 
thermocouple of unequal diameters that does not require compensation at lower 
fluid temperature frequencies nor any knowledge of the fluid velocity or 
properties. 4 - 5 The results of experimental measurements are presented along with 
the suggested procedure for the reduction of the data from the three wire 
thermocouple as shown in fig. 1. 
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2. THEORY 


Use of a single wire thermocouple requires knowledge of the fluid velocity 
and properties to determine the natural frequency. The latter quantity is necessary 
to establish a frequency dependent compensation factor for the measured signal. 


2.1 Natural Frequency of Single Wire 

The conservation of energy for a single wire is : 


^ = <o„(T,-T) 


( 1 ) 


where the natural frequency 


Here: 


(On 

T 

Tg 
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h 
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CO. 


4h 

pcD 


( 2 ) 


= natural frequency (sec* 1 ) 

= temperature of the thermocouple (°K) 

= temperature of the gas (°K) 

= time (sec) 

= wire density (kg/m 3 ) 

= heat capacity of wire (J/kg - *K) 

= heat transfer coefficient (W/m 2 - 'K) 

= wire diameter (m) 

It is convenient to rewrite the heat transfer coefficient in terms of the Nusselt 
number 


h = 


— Nil 
D 


(3) 
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where 

kg = thermal conductivity of the gas (W/m - *K) 

Nu = Nusselt number. 


Thus, the natural frequency in Eq. (2) becomes 


©n 


4k g Nu 
(pc)D 2 ' 


(4) 


where 6 

Nu = 

C c Re m Pr 1/3 - 

(5) 

and 

Re = 

VD 



Pr = 




Here, Re and Pr are the Reynolds and Prandtl numbers, respectively and 
fi = viscosity of gas (kg/m-sec) 
t)g = Kinematic viscosity of gas (m 2 - sec 1 ). 

V = velocity of gas (m/s) 


The natural frequency of a single wire from Eqs. (4) and (5) can now be written 
in the form. 




k g C c 4Re m Pr 1/3 
(pc)D 2 


or separating the wire diameter 


( 6 ) 
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(7) 


co n = CD™ -2 

where the coefficient C depends on both fluid and wire parameters or 


4k CV m Pr 1/3 
C = — 8 0 


(pc)t) 


(8) 


Substituting values corresponding to a type K thermocouple and a gas velocity V = 
25 m/s for air at standard conditions into Eqs. (7) and ( 8 ), one obtains values of the 
natural frequency: 

co nfsec 1 ) D (m) D(mil) 

52.6 50.8x1 O ' 6 2 

18.3 101. 6x1 O ' 6 4 

6.5 203.2x1 O ' 6 8 


2.2 Determination of Gas Temperature 

The use of a single wire thermocouple requires knowledge of the fluid 
velocity and additional properties such as viscosity or density. There are three 
unknowns that appear in Eqs. (1) and (7): a) the coefficient C given by Eq. ( 8 ) b) m the 
exponent of the wire Reynolds number in Eqs. ( 6 ), (7) and ( 8 ) and c) the gas 
temperature T g . 

Consider three thermocouples of unequal diameters Dj < D 2 < D 3 as shown in 
fig. 1. The conservation of energy equation for each wire reduces to 


dT, 

dt 


CD ™- 2 (Tg-T,) 


(9) 


dt 


CD ™" 2 (T g -T 2 )- 


( 10 ) 
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( 11 ) 


dt 


CDr 2 (T -T 3 ). 


In principle it is possible to solve for the three unknowns that appear in Eqs. 
(9)-(ll) in terms of the instantaneous values of the wire temperatures Ti, T 2 and T 3 
and their derivatives. However, because of the empirical nature of the Nusselt 
number given in Eq. (5) and uncertainty in the wire diameters near the 
thermocouple junctions, we choose to eliminate the coefficient C and solve for the 
gas temperature Tg from Eqs. (9)-(ll). Moreover, because the conservation of energy 
equations are linear at constant ambient velocities, we choose to use transformed 
data profiles. 

Thus, taking the Fast Fourier Transform (FFT) of Eqs. (9)-(ll), squaring Eq. (10) 
and dividing by the product of Eqs (9) and (11), one obtains 


T 2 ') 

i 2 


f Dl ) 

(T g — T 2) 2 

< T ^ 3 ) 



(Tg-T.KT.-T,) 


( 12 ) 


Substituting r = T 2 2 / T 3 T, 


(13) 


into Eq. (12) and anticipating that the wire diameter ratio D 2 70^3 
or 


2 A 


D 




m-2 


= 1+e, 


the transform of the gas temperature becomes 


0 ( 1 ) 


(14) 


— -b - Jb 2 -4a c 

rp O V 0 O O 

8 " 2a] 


(15) 
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Here, 


b D — b + 2 g T 2 +0(g) 2 
-4a 0 c 0 = 4 e T 2 2 (r-l)+0(e) 2 
b o _4a o c o = b 2 +Gy+0(e) 2 
and a o = r-l-G 

where b = 2T 2 — r(T 3 +T,) (16) 

and y — 4T 2 (b+T 2 (r-l)). 


We now seek a Taylor series expansion of Eq. (15) for the transformed gas 
temperature about the point e = 0. Thus, to first order in e where 



f(s) = -4a ri c o = ^/b 2 + e y 


and 

, . 1 1 
g(e) = = , 

a 0 r-l-e 


one obtains 

f(e) = b + — g+0(g) 2 
2b 

(17) 

and 

g(e) = l+4+0(e) 2 
a a 

(18) 

where 

a = r-1. 

(19) 


Substituting Eqs. (17) and (18) into Eq. (15) and retaining terms O(e), one obtains 


b 2T 2 T 2 b^l 
--e — + + — 

ci \ ^ J 


where a, b are defined by Eqs. (13), (19) and (16), respectively. 


( 20 ) 
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We now define the transform of the gas temperature for the particular case 
where the wire diameter ratio D* /DjD., = 1 in the form 


T g (0) = - 


b 

a 


( 21 ) 


or 


T„(0) 


T ; (T,T,_+ T 2 T 3 _- 2T,T 3 ) 
(T; - T,T,) 


( 22 ) 


Substituting for the quantity b defined by Eq. (21) into Eq. (20), one obtains the final 
result 


T g (e) = T g (0) + e 


‘ X ' 
v aT 8 (0) ; 


(t 2 -t 8 (o)) j 


(23) 


V 

where a = — — — 1 

TT 

A 1 x 3 

and T g (0) is defined by Eq. (22). The inverse transform of T g (e) becomes the 
reconstructed gas temperature or 


T g = FFT-'[T g (e)] (24) 

for assumed values of e . 

It is useful to consider the effects of wire diameter on the magnitude of the 
perturbation parameter e . In general it is desirable to have small diameter wires so 
that their natural frequencies are shifted to larger values. Listed in table 1 are three 
possible wire combinations that are available commercially along with the 
magnitude of e . The values of e are computed from Eq. (14) assuming that the wire 
diameters are exactly as listed by the manufacturer which is very unlikely in 
practice. For example, a typical value of the exponent m = 1/2 in Eq. (14) suggests 
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that a 7% error in D2 would yield a value of e < 0.2 for the wire diameters listed in 
the first row of Table 1. 

TABLE 1 - Wire Diameter Combinations (m = 1/2) 

Di (mil) P?(mil) Ch (mil) § 

2 4 8 0.0 

2 3 5 0.17 

2 3 4 - 0.16 

3. EXPERIMENT 
3.1 Apparatus 

In the present experiment, thermocouple sensors are exposed to a constant 
velocity air stream (<25 m/s) of varying temperature. In particular, the dynamic 
response of the thermocouple is measured for a periodic temperature profile of 
varying frequency. A rotating wheel configuration is used to deliver the test air 
stream to the proposed sensors. A similar experimental apparatus was described in 

detail by Elmore et al 7 and Forney et al. 8 

A schematic of the rotating wheel apparatus used in the present experiment is 
shown in Fig. 2. As the wheel rotates, holes pass the two air supply tubes (3/4-in. i.d. 
copper) that allow slugs of hot (~55*C) and cold (~30*C) air to alternately enter a 
transition tube assembly mounted directly above the rotating wheel. In the 
transition tube the slugs of hot and cold air coalesce into a single air stream 
providing a periodic temperature profile covering a range of gas temperature 
frequencies from roughly 1 to 60 Hz. A 3/4-1 /2 in. smooth copper adapter was 
inserted at the end of the coalescing tube (location of thermocouple) in Fig. 2 to 
increase the air velocity to 25 m/s. 
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The analog temperature signal is digitized with a Data Translation DT-2801 
A/D board mounted in an expansion slot of an IBM AT compatible computer as 
shown in Fig. 2. ASYST software was loaded onto the hard disk of the personal 
computer and this provided a flexible system for data storage, manipulation, and 
display. The true temperature profile of the airstream is measured with a constant 
current anemometer (TSI 1054- A) and sensor (1210-T1.5). 

The three type K thermocouple wires were of diameter D} = 50.8pm (2 mil), 
E >2 = 101.6pm and D 3 = 203.2pm as shown in fig. 1 . The thermocouple junctions 
were fabricated by the Research Instrumentation Branch at the NASA Lewis Center. 
The wires were cut with a razor blade to produce a flat edge perpendicular to the 
wire axis. The wire segments were then mounted in a fixture such that the two 
faces of the junction were held together by springs. The junction was then formed 
by laser heating. The laser used was a KORAD model KWD Nd: YAG laser, 
operating at a wave length of 1060 nm. The power settings depend on the wire size 
and material, but for the thermocouples described in this report, a pulse duration of 
4 ms was used delivering a total energy of approximately 2J. 

3.2 Procedure 

ASYST software was developed to acquire temperature data sequentially from 
the thermocouples and constant current anemometer. The data are digitized for 
four channels at a sampling rate of 512 Hz per channel for a total sample time of 0.5 
s. After data acquisition, the Fast Fourier Transform (FFT) is taken for each channel. 
The largest peak (first harmonic) in the amplitude of the FFT is located for each 
channel which provides the frequency of the temperature profile. The ASYST 
software also records the amplitude ratio of the first harmonic for the 50.8pm 
thermocouple - to - anemometer output. The amplitude ratio was recorded for each 
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frequency setting of the rotating wheel and provided a value for the natural 
frequency o>i of the small wire. 

The FFT of the signal from each thermocouple channel was substituted into 
Eq. (19) by the ASYST code along with an estimate of the parameter e . The inverse 
transform of T g by the ASYST code provided a reconstruction of the true gas 
temperature. 

4. RESULTS AND DISCUSSION 

The temperature profiles for the ambient gas (large amplitude) and the three 
type K thermocouple wires are shown in fig. 3. The three wire diameters used in the 
present study are Dj = 50.8pm (2mil), D 2 = 101 . 6 pm and D 3 = 203.2pm. The angular 
frequency of the gas temperature in fig. 3 is co = 86 s ' 1 (~ 14 Hz) while the natural 
frequency of the small wire coi = 42 s _1 . The natural frequency of the small wire 
corresponds to an ambient gas velocity of roughly 10 m/s. 

Figures 4, 5 and 6 illustrate the effects of increasing the magnitude of the wire 
diameter parameter e defined by Eq. (14). The small phase angle between the large 
ambient gas temperature profiles that appear in figs. 4 and 5 and the smaller 
temperature profiles reconstructed from the thermocouple probe allow one to 
distinguish between both profiles in fig. 6 . In the latter case, the temperature profile 
on the left is that of the ambient gas while the profile on the right was determined 
by taking the inverse FFT of T g (e) defined by Eq. (23) with e = 0.2. In all three 

figures the natural frequency of the small wire coi = 61 s * 1 (V ^ 20 ms' 1 ) while the 
angular frequency of the ambient gas temperature is co = 209 s _1 (~ 33 Hz). 

It is clear that a value of the parameter e = 0.2 provides a signal that requires 
no compensation for the frequency ratio co/coi = 3.4. It should also be noted that a 7% 
error in the diameter of D 2 would account for a value of e = 0.2. Further, the 
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coefficient of e in Eq. (23) depends on the difference between the uncompensated 
FFT of the signal from the wire of diameter D 2 and the compensated FFT of the 
signal from the three wire probe with a particular diameter ratio D* = D, D 3 (e = 0). 
Thus, the contribution to the signal by the second term on the right of Eq. (23) is 
negligible at low frequencies but increases with larger values of co/coi. 

The amplitude ratio of the three wire thermocouple output-to-true gas 
temperature as a function of gas temperature frequency has been plotted in fig. 7. 
Also shown is the first order response of the small wire. The open symbols 
correspond to a value of e = 0 in Eq. (23) where the amplitude ratio requires no 
compensation for go < 2coi. For the range of values 2coi < go < 5Grtbe amplitude ratio 
follows the function 


A 


1.92 


( y 103 

1 (0_ ' 

) 


(25) 


It should be noted that the first order response is also proportional to (go/goj) - 1 at 
large frequencies. 

A significant improvement is observed, however, when the parameter e = 0.2 
as indicated by the solid symbols in fig. 7. For the latter case, no compensation is 
required for (o<5coi. As discussed earlier, the amplitude ratio data coelesce at low 
gas temperature frequencies where the coefficient of e in Eq. (23) approaches zero. 

Both values of e = 0.0 or 0.2 are of little value, however, when go >6goi where 
the amplitude ratio for the three wire thermocouple reduces to the first order 
response of the smaller wire. The latter result appears to occur at a gas frequency of 
~ 60 Hz. For the present experimental system the signal - to - noise ratio may be 
responsible for the small response at large frequencies. This situation may be 
improved by replacing the two larger 
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wires of 4 and 8 mil diameter by wires of 3 and 4 mils, respectively, as listed in table 

1 . 

5. CONCLUSIONS 

Theory and experimental measurements are compared with a novel three 
wire thermocouple. Signals from three wires of unequal diameters are recorded 
from the thermocouple suspended in constant flow with a periodic temperature 
fluctuation. It is demonstrated that the reconstructed signal from the three wire 
thermocouple requires no compensation for 0)<5(0i where oo i is the natural 
frequency of the smaller wire. The latter result represents a significant 
improvement compared to previous work with two wire thermocouples. A 
correction factor has also been derived to account for wires of arbitrary diameter. 
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Fig. 2 - Schematic of rotating wheel apparatus. 
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Fig. 3 - Temperature profiles for the gas (large amplitude) and three type 
K thermocouple wires. Three wires: Di = 50 . 811 m, toi = 42 s' 1 , D 2 : 
101.6pm, D 3 = 203.2pm. Gas temperature frequency: co = 86 s' 1 . 
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Fig. 4 - Gas temperature profile compared to reconstructed temperature 
derived from three wire thermocouple. ol>i = 61 S' 1 , to = 209 s* 1 . 
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Fig. 6 - 
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Fig. 7 - Amplitude ratio of three wire thermocouple output- to- true gas 
temperature. Also shown is first order response of small wire. 


III. MULTIWIRE THERMOCOUPLES IN 
REVERSING FLOW 
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1. INTRODUCTION 

The advantages of thermocouples are their low cost, reliability and simplicity 
since they do not require optical access or elaborate support electronics. The design 
of a thermocouple represents a compromise between accuracy, ruggedness and 
rapidity of response. Although thermocouples are suitable for the measurement of 
high frequency temperature fluctuations (< 1 KHz) in a flowing gas or liquid, the 
measured signal must be compensated since the frequency of the time-dependent 
fluid temperature can be much higher than the natural frequency of the 
thermocouple probe. 1 Moreover, use of a single wire thermocouple requires 
knowledge of the fluid velocity (normally assumed constant) and fluid properties 
(e.g., viscosity, density, etc.) to determine the convective heat transfer coefficient of 
the wire and its natural frequency. 

The present paper describes the performance of novel two or three wire 
thermocouples with unequal wire diameters for use in unsteady or irregular flows 
of varying velocity. In this case the time constant or natural frequency of the 
thermocouple wire is time dependent and the normal procedure for compensation 
of high frequency signals is impractical. Previous work 2 3 ' 4 indicated the usefulness of 
the concept in constant velocity flows where it was demonstrated that no 
compensation was required nor any knowledge of the fluid velocity or properties 
over a useful range of fluid temperature frequencies. In the present paper the 
multiwire thermocouple has been tested in a reversing flow field. The results of 
experimental measurements are presented along with the suggested procedure for 
the reduction of the data from the multiwire thermocouple as shown in fig. 1. 


2. THEORY 

Unsteady fluid motion develops in a variety of circumstances. For example, 

surfaces of discontinuity behind sharp edges such as airfoils and adverse pressure 
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gradients on blunt objects and in sharp corners lead to flow separation and the 
formation of recirculating eddies. 5 Other examples of irregular fluid motion are a.) 
cellular vorticies in stratified fluids or b.) oscillating flows in the exhaust from 
certain engines (e.g., Stirling, etc.). 


2.1 Natural Frequency of Single Wire 

The conservation of energy for a single wire is 


dT 


= (0„ (T. - T) 


dt ' 8 

where the natural frequency 


( 1 ) 


Here: 


(0 n = 


4h 

pcD 


( 2 ) 


(On = natural frequency (sec 1 ) 

T = temperature of the thermocouple OK) 

Tg = temperature of the gas OK) 

t = time (sec) 

p = wire density (kg/m 3 ) 

c = heat capacity of wire (J/kg - *K) 

h = heat transfer coefficient (W/m 2 - °K) 

D = wire diameter (m) 

It is convenient to rewrite the heat transfer coefficient in terms of the Nusselt 
number 


h = —Nu 
D 


( 3 ) 


where 


k g 


thermal conductivity of the gas (W/m - °K) 
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Nu = Nusselt number. 
Thus, the natural frequency in Eq. (2) becomes 


CO 


n 


4h g Nu 
(pc)D 2 ' 


where* Nu = C D Re m Pr 1/3 - 

Re = M 


(4) 

(5) 



Here, Re and Pr are the Reynolds and Prandtl numbers, respectively and 
|i = viscosity of gas (kg/m-sec) 
u g = Kinematic viscosity of gas (m 2 - sec- 1 ). 

V(t) = time dependent fluid velocity (m/s) 

The natural frequency of a single wire from Eqs. (4) and (5) can now be written 
in the form. 


g o 

co n = -S 


k„ CL 4Re m Pr 1/3 


(pc)D 2 

or separating the wire diameter and fluid velocity 


to n = C|V| m D m “ 2 

where the coefficient C depends on both fluid and wire parameters or 


(6) 


(7) 
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c = 


( 8 ) 


4k g C p Pr 1/3 
(pc)u™ 

Substituting values corresponding to a type K thermocouple and a maximum gas 
velocity V max = 25 m/sec 1 for air at standard conditions into Eqs. (7) and (8), one 
obtains values of the maximum natural frequency G> max : 
fflmay(sec? l) D (m) P(nu ll 

52.6 50.8x10-6 2 

18.3 101.6x10-6 4 

6.5 203.2x10-6 8 


2.2. Two Wires 

The conservation of energy given by Eq. (1) for a single wire contains two 
unknowns: co n the natural frequency and the desired gas temperature T g . Each 
additional wire of a different diameter adds an unknown natural frequency co n . We 
also assume for a multiwire thermocouple that additional information is available 
concerning the ratios of natural frequencies. The latter provides a sufficient number 
of independent equations to determine the time-dependent gas temperature T g . 

Consider two thermocouple wires of unequal diameter Dj < D 2 . The 
conservation of energy for each wire becomes 


dT\ 

dt 


C|V| m D5 n-2 (T g — Tj) 


(9) 


dT 2 

dt 


C|V| m D?" 2 (T g -T 2 ). 


( 10 ) 


In principle, knowledge of the wire diameters provides an additional equation 
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m-2 


(ID 


= fiV 

(0 2 V^2> 

representing the ratio of natural frequencies as defined by Eq. (7). Dividing Eq. (9) by 
Eq. (11) and solving for the gas temperature one obtains 



where Tj = dTj/ dt. 

Our experience with multiwire thermocouples in constant flow 3 - 4 and the 
present study indicate that construction of the gas temperature from Eq. (12), which 
requires the measured wire signals Ti, T 2 and their derivatives, is unstable and 
unreliable. We have found in all cases, however, that satisfactory results are 
obtained by taking the Fourier transform of the conservation equations and 
deriving the gas temperature in terms of transformed variables. For the case of 
unsteady flow, independent knowledge of the gas velocity is therefore required. 

Thus, taking the Fast Fourier Transform (FFT) of Eqs. (9) and (10) and 
dividing the expressions, the transform of the product of gas temperature and fluid 
velocity is obtained in the form 
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where Ti v = FFT(|V| m Ti) and T gv = FFT (|V| m T g ). For convenience we also define 
the parameter (J = coi /(02 or 


P = 


^Di ^" 2 




(14) 


The inverse transform of T gv provides the reconstructed gas temperature or 



2.3. Three Wires 

Consider three thermocouple wires of unequal diameters Di < D 2 < D 3 as 
shown in fig. 1 . The conservation expression for each wire becomes 


dT 


dt 




(15) 


dT 2 

dt 


= ® 2 (T g -T 2 ) 


(16) 


dlj 

dt 


- «3( T g - T 3) 


(17) 


with the additional constraint 


Cjl) 7 

a = — - 


01CO3 


< p| ^ m_2 


(18) 


Taking the FFT of Eqs. (15)-(17) / squaring Eq. (16) and dividing by the product 
of Eqs. (15) and (17), one obtains 
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( 19 ) 


f "j _ a ( T gv T 2v ) 

^ T1T3 J ^Tg V — T| v )(Tg V - T 3 v j 

Substituting r = T 2 2 /T 3 T 1 and solving for T^, one obtains 


-b-> 2 * * * -4ac qo) 

** 2a 

where a = r - a 

b = 2aT 2v -r(T lv +T 3v ) 
c = rT lv T 3v — aT 2v 

and = FFT(|V| m T g ), T lv = FFT^|V| m T 1 j. The inverse transform of T gv provides 
the reconstructed gas temperature or 



( 21 ) 


3. EXPERIMENT 

3.1 Apparatus 

In the present experiment, thermocouple sensors are exposed to a reversing 

velocity air stream (<25 m/s) of varying temperature. In particular, the dynamic 

response of the thermocouple is measured for a periodic temperature profile of 
varying frequency. A rotating wheel configuration is used to deliver the test air 

stream to the proposed sensors. A similar experimental apparatus was described in 
detail by Elmore et al 7 and Forney et al 8 . 
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A schematic of the rotating wheel apparatus used in the present experiment is 
shown in Fig. 2. As the wheel rotates, holes pass the two air supply tubes (3/4-in. i.d. 
copper) that allow slugs of hot (~55°C) and cold (~30°C) air to alternately enter the 
tube assembly mounted directly above the rotating wheel. The slugs of hot and cold 
air form a reversing air stream providing a periodic temperature profile covering a 
range of gas temperature frequencies from roughly 1 to 60 Hz. 

The analog temperature signal is digitized with a Data Translation DT - 2801 
A/D board mounted in an expansion slot of an IBM AT compatible computer as 
shown in Fig. 2. ASYST software was loaded onto the hard disk of the personal 
computer and this provided a flexible system for data storage, manipulation, and 
display. The true temperature profile of the airstream is measured with a constant 
current anemometer (TSI 1054- A) and sensor (1210-T1.5). 

The three type K thermocouple wires were of diameter Di = 50.8pm ( 2 mil), D 2 
= 101.6pm and D 3 = 203.2pm as shown in fig. 1. The thermocouple junctions were 
fabricated by the Research Instrumentation Branch at the NASA Lewis Center. The 
wires were cut with a razor blade to produce a flat edge perpendicular to the wire 
axis. The wire segments were then mounted in a fixture such that the two faces of 
the junction were held together by springs. The junction was then formed by laser 
heating. The laser used was a KORAD model KWD Nd: YAG laser, operating at a 
wave length of 1060 nm. The power settings depend on the wire size and material, 
but for the thermocouple described in this report, a pulse duration of 4 mn was used 
delivering a total energy of approximately 2J. 

3.2 Procedure 

ASYST software was developed to acquire temperature data sequentially from 
the thermocouple and constant current anemometer. The data are digitized for four 
channels at a sampling rate of 512 Hz per channel for a total sample time of 0.5 s. 
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After data acquisition, the Fast Fourier Transform (FFT) is taken for each channel. 
The inverse transform of T gv by the ASYST code provided a reconstruction of the 
true gas temperature. 

4. RESULTS AND DISCUSSION 

The velocity profile of the airstream was measured with the constant current 
anemometer. Figure 3 illustrates the velocity profile for a wheel frequency of 10 Hz. 
For the positive velocities in fig. 3 the gas reached a temperature of roughly 25 to 
30°C corresponding to the cold slugs while the negative velocity components were 
hot air slugs of approximately 50°C. 

In the present study the gas velocity and temperature were not measured 
simultaneously. Rather, it was assumed that the temperature measured with the 
anemometer was proportional to the gas velocity in the reversing flow field or V(t) 
a Tg(t). A typical data profile is illustrated in fig. 4. The uniform profile on the right 
is the temperature profile measured with the anemometer. The temperature profile 
on the left was reconstructed with Eg. (20) assuming that the exponent m = 0.5 in 
Eqs. (7) and (21) and a value of the diameter ratio a = 1.0 for the three wire 
thermocouple. 

The value of the exponent m = 0.5 was chosen for the Reynolds number that 
appears in the Nusselt number correlation Eq. (5) and the wire natural frequency wn 
given by Eq. (7). At the maximum gas velocity of 25 m/s the maximum Reynolds 
number for the three wires in the thermocouple assembly are 80, 160, 320 
corresponding to the smallest to largest diameter wire, respectively. Since the 
recommended values of the exponent m are 0.5 for Re > 40 and 0.4 for 1 < Re < 40, 6 
we chose a constant value m = 0.5 in the present study. 
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4.1 Two Wires 

The amplitude ratio of thermocouple-to-gas temperature was determined 
using a two-wire combination with Di= 50.8 pm and D2 = 101.6 pm from the 
schematic in fig. 1. The gas temperature was reconstructed with Eq. (13) assuming m 
= 0.5. The amplitude ratio was determined by averaging peak temperatures over six 
or eight cycles. The results are indicated in fig. 5 as a function of the parameter P = 
©1 / ©2 = (Di /D 2 )- 3 / 2 . In fig. (5) the natural frequency ©1 = 59 s- 1 is defined as the 
frequency at the maximum gas velocity for the small wire Di = 50.8 pm. It is clear 
that by amplitude ratio is strongly affected by p and that the optimum value is P = 
2.55. 

Figure 6 illustrates amplitude attenuation as a function of gas temperature 
frequency for fixed P = 2.55. The two wire combination requires no compensation 
for frequencies ©/ ©1 < 2.3 corresponding to frequencies < 22 Hz in the reversing 
flow field. 

4.2 Three Wires 

The amplitude ratio of thermocouple-to-gas temperature was determined for 
the three wire combination as shown in the schematic of fig. 1. The gas temperature 
was reconstructed with Eqs. (20) and (21) with the exponent m = 0.5. The results are 
indicated in fig. 7 as a function of the parameter a=(©2/©,© 3 )=(D2/D 1 D 3 )“ 3/2 . It is 
clear that the amplitude is strongly affected by a and that the optimum value is a = 
1 . 0 . 

Figure 8 illustrates amplitude attenuation as a function of gas temperature 
frequency for fixed a = 1.0. The three wire combination requires no compensation 
for frequencies ©/ ©1 < 3.5 corresponding to frequencies < 34 Hz in the reversing 
flow field. 
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5. CONCLUSIONS 


Measurements are recorded for multiwire thermocouples consisting of either 
two or three wires of unequal diameters. Signals from the multiwire probe are 
recorded for a reversing gas flow with a periodic temperature fluctuation. It is 
demonstrated that the reconstructed signal from the multiwire thermocouple 
requires no compensation provided o )/ (Oi < 2.3 for two wires or G)/(fli < 3.6 for three 
wires where coi is the natural frequency of the smaller wire based on the maximum 
gas velocity. The latter results were possible provided Fourier transformed data 
were used and knowledge of the gas velocity is available. 


61 



6. REFERENCES 


1. M. D. Scadron, and I. Warshawski, Experimental Determination of Time 
Constants and Nusselt Numbers for Bare-Wire Thermocouple in High 
Velocity Air Streams and Analytic Approximation of Conduction and 
Radiation Errors, NACA TN-2599 (1952). 

2. G. C. Fralick, Patent Disclosure, April, 25 (1990). 

3. Forney, L. J and Fralick, G. C., Two Wire Thermocouple: Frequency Response 
in Constant Flow, Rev. Scientific Instruments 65, 3252 (1994). 

4. Forney, L. J. and Fralick, G. C., Three Wire Thermocouple: Frequency 
Response in Constant Flow, Rev. Scientific Instruments (submitted, 1994). 

5. Prandl, L., Essentials of Fluid Dynamics, Hafner Publishing Co. (1952). 

6. Incropera, F. P. and D. P DeWitt, Fundamentals of Heat and Mass Transfer 
2nd ed. New York: John Wiley and Sons, 1985. 

7. D. L. Elmore, W. W. Robinson, and W. B. Watkins, Further Development of 
the Dynamic Gas Temperature Measurement System: Vol. I-Technical Efforts, 
NASA CR-179513 (1986). 

8. L. J. Forney, E. L. Meeks, J.Ma and G. C. Fralick, Measurement of Frequency 
Response in Short Thermocouple Wires, Rev. Sci. Instrum 64 3236 (1993). 


62 


7. FIGURES 



Schematic of three wire thermocouple 
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Figure 2 






03 x 



66 










oiiva Bamiidwv 


71 


Figure 8 - Amplitude ratio of three thermocouple output versus gas 

temperature frequency, toi = 59 s' 1 . 



IV. APPENDIX - ASYST CODES 
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1. ASYST CODE -CONSTANT FLOW 


\ FF3P.F0K 
ECHO . OFF 
FORGET . ALL 

REAL DIM[ 512 ] ARRAY TIKE 

COKF'LEX DIKE 512 ] ARRAY TEAKS 1 

COMPLEX DIME 512 ] ARRAY TRAN 32 

COMPLEX DIME 512 ] ARRAY TRANS3 

REAL DIME 10 ] ARRAY ANS 

REAL DIME 512 ] ARRAY ZMAGO 

REAL DIME 512 ] ARRAY 2ARG0 

INTEGER SCALAR NUM 

REAL SCALAR MEN 

REAL SCALAR KEN1 

REAL SCALAR MENS 

REAL SCALAR MENS 

REAL SCALAR MEN 4 

REAL SCALAR CD- 

REAL SCALAR TI 

REAL DIME 512 ] ARRAY FREQS 

REAL DIME 512 ] ARRAY SIGNAL 

INTEGER DIME 512 , 4 ] ARRAY DATA. BUFFER 


25 CD : - 

. C C DC DATA . BUFFER : = 

LOAD . OVERLAY ACQU I S . SOV 
0 3 A/D. TEMPLATE DEMO .TEMPLATE 

DATA . BUFFER TEMPLATE . 3UFFER 
CYCLIC 

512. TEMPLATE . REPEAT 
CD CONVERSION. DELAY 
DEMO . TEMP LATE A/D . I N I T 

\ SYNCHROS' I ZE 
A/D. IN --ARRAY 

LOAD . OVERLAY’ WAVEOFS . SOV 

512 REAL RAM? 1-4. * CD t 1000. 



DATA. BUFFER XSECTE 1 ] 2048 - 
Signal : — i on a. _ 

—.9336 ** 1324. * 

SIGNAL := SIGNAL 
.05 SET . CUTOFF . FREQ 
SMOOTH SIGNAL : - SIGNAL 
MEAN MEN := SIGNAL 
KEN - SIGNAL := 

TIME SUE! 100 , 300 , 1 ] SIGNAL 
XY . AUTO . PLOT 

SIGNAL 



\ CENTER ON ORIGIN 
SUE[ 100 . 300 , 1 ] 
\ PLOT CHANNEL 0 
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FFT \ TAKE FFX 

TKANS1 := TRAN’S 1 ZMAG ZMAGO : = \ MAGNITUDE 

5 SET.*. PCI NTS 
1 SET.#. OPTIMA 

ZMAGO SUB [ 1 , 120 , 1 ] LOCAL. MAXIMA \ 

SWAP NUM := NUM \ 

1 CD / 4 . / 512. / 1000. * PI * 2. * \ 


ANS 

[ 1 J := ANS I 1 ] 

DROP ANS [ 2 3 

DEG 

TRAN 31 ZARG ZARGO : 

= ZARGO [ NUM ] 

ANS 

[33: = 

TRAN 

El [ NUM ] DEG ZARG 

ANS [ 3 ] := 


OF FFT 


FIND INDEX AND MAX OF MASK 
INDEX OF MAX MAGNITUDE 
FREQUENCY AT MAX MAGNITUDE 

FIND ARGUMENT OF FFT 
ARGUMENT AT MAX MAGNITUDE 


DATA. BUFFER XSECT[ 2 ] 204S 
SIGNAL : = SIGNAL 
. 04725 * 


\ CHANNEL 1 ON STACK (THERMOCOUPLE ; 


\ SMOOTH SIGNAL 
MEAN MEN 1 := 


ilGNAI 

SIGNAL MEN1 - SIGNAL := 


\ SMOOTH DATA ( CYCLES/PC INT ) 


UB: _ 100 . 300 , 1 ] SIGNAL SU3[ 100 , 300 , 1 ] 


XY . DATA . PLOT 

SIGNAL 
Fr * 

TRAN 31 : - TRAN SI ZMAG ZMAGO : = 

ZMAGO SUE! 1 , 120 , 1 ] LOCAL. MAXIMA 
SWAP NUM : = NUM 

- ~ CL’ / 4. / 512. / 1000. + FI ♦ 2. * 
AN2 C < T : = ANS [ 4 ] DROP ANS [ 5 ] ; 

- ^*A X w . iARG .. Ar.CU : - ZARGO [ NUM 3 
ANS : S : := 

TRANS 1 [ NUM 3 DEG ZARG ANS [ 0 ] := 


\ PLOT CHANNEL 1 


\ TAKE FFT 
\ MAGNITUDE OF FFT 

\ FIND INDEX AND MAX OF MAC 
\ INDEX OF MAX MAGNITUDE 
\ FREQUENCY AT MAX MAGNIFY! 

\ ARGUMENT AT MAX MAIN I TUTS 


DATA. BUFFER XSSC 
SIGNAL := SIGNAL 
.05457 * 

SIGNAL := SIGNAL 


] 2C4S 


\ CHANNEL 2 ON STACK (THEE 


\ &MOOTH SIGNAL := SIGNAL \ SMO^iTH DA~A (CY'“" T ' = 'S'" i 

MEAN MEN 2 := SIGNAL MEN2 - SIGNAL := \ CENTER 0N*0RIGIn“"" 

300 . 1 ] SIGNAL SUB[ 100 , 300 , 1 ] 

\ PLOT CHANNEL 2 


TIME sub: ioo 
XY . DATA . PLOT 


SIGNAL 

FFT 

THANS: 


\ TRANSFORM OF CHANNEL 
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\ CHANNEL 3 ON STACK (THERMOCOUPLE) 


DATA. BUFFER XSECT[ 4 ] 2048 - 
SIGNAL : = SIGNAL 
.04725 * 

SIGNAL := SIGNAL 

\ SMOOTH SIGNAL := SIGNAL \ CEN~^F ON^OR'G'N^^' 1 ' 0 *’' * " 

MEAN KENS := SIGNAL MENS - SIGNAL := \ CEN.nR ON CR.C.N 

TIKE SUB[ 100 , 300 , 1 ] SIGNAL SUBl 100 , 30^ 

XV. DATA. PLOT 


SIGNAL 

FFT 

TRANSS 


\ PLOT CHANNEL 3 
\ TRANSFORM OF CHAN N El 


TRANS 2 TRAN Si * TRANSS TRANSS * + 

2 . TRAN Si * TRAN S3 * - TRANSS * 
TRANS 2 TRANS 2 * TRAN31 TRAN S3 * - / 
TRAN Si : - 


TRAN Si 2 MAG 2MAG0 := 


\ 

ZMA 

vjO iJyJZ 

I i 

, 120 

\ 

ANS 

: 7 ] 

• = 


TRAN Si 




2REAL S 

IGNAL 

: = 

S i. Gi'A 


MEAN 

MEN 4 : 

r c , 

IGNAL 


SMOOT 

K SIGN 

AL 

: = 

C 7- 

O J. 

GNAL 

SIGN 

AL 

: - 

XI 

yj'p c;? 

El 100 


SCO , 


GNAL 

sub: ioo 

, SCO 

XY 

' . DATA 

T)T 



CR 




\ 


FREQ 


MAG 


ANS [ 

1 ] 


CR 


ANS [ 

5 ] 

ANS 

1 L — . 

V 

\ 

ANS 

[ 7 ] 

2 

ANS [ 

\ 

ANS 

[ 7 ] 

AN 2 

i [ 2 : 

ANS 

[ 3 j 

■i 

— X 

* ANS 


ANS 

r *3 "i 

L ^ J 

_ 1 

* ANS 


MEN 

? MEN j 

1 ? 

MENS # 


LOCAL. MAXIMA 


PHI 


CR 


\ CON STRUCT S a. GN Ai. 


\ AMPLITUDE CF FIRST HA; 


INVERSE 


r r* 


\ 


\ ANGULAR FF.EG CHANNEL 0 
\ AMPLITUDE RATIO CKANNEI 


GNA- 


! /O 


CR 


CR 


\ FORGET. ALL 
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2. ASYST <X)DE - REVERSING FLOW 


\ REV. FOR 
ECHO . OFF 
FORGET. ALL 


REAL DIME 

512 

1 

J 

ARRAY 

TIME 

COMPLEX DIME 

51 

2 ] ARRAY TRAN 

REAL DIM' 

512 

] 

ARRAY 

SIG1 

REAL DIM! 

512 

] 

ARRAY 

SIG2 

REAL DIM[ 

512 


ARRAY 

SIG3 

REAL DIM[ 

512 

J 

ARRAY 

DIF1 

REAL DIM: 

512 

■j 

ARRAY 

DIF2 

REAL DIM: 

512 

] 

ARRAY 

DIF3 

REAL DIM: 

512 

J 

ARRAY 

R1 

REAL DIM' 

1C 

] ARRAY AX 3 

REAL DIM: 

512 

J 

ARRAY 

ZMAGO 

REAL DIM: 

o 

] 

ARRAY 

LARGO 


INTEGER SCALAR N’JM 
REAL SCALAR MEN 
RiAL SCALAR KEN1 
REAL SCALAR MENS 
REAL SCALAR MENS 
REAL SCALAR KEN 4 
REAL SCALAR CD 
KEAi_: SCALAR AL 

HEAL DIM-: 512 j ARRAY SIGNAL 

i N _ Ffir n T.~ 'S' F. * • ~ A TE A V 


or . _ 

.OCC: DATA . ELTFER 
LOAD. OVERLAY ACQ’JII 
0 3 A/D . TEMPLATE 




- fi-. n . . 


'FEE TEMF LATE . BUFFER 


£12. TEMPLATE . REPEAT 
CL CONVERT I ON . DELAY 
DIMS . TEMPLATE A/L . I N I T 
\ i. SYNC. PERIOD 
\ SYNCHRONIZE 
A/L. I N: ARRAY 

LOAD. OVERLAY WAVEOPS.SOV 

512 REAL RAM? 1-4. + CD * 1000. / TIME := 


\ DEFINE AN A/L TEMPLATE 
\ DECLARE ARRAY AS A TEMFLATE BUFFER 
\ SET TEMPLATE BUFFER TO CYCLIC MODE 
\ SET REPETITIONS FOR I/O INSTRU^Lo'.' 
\ CONVERSION RATE ( MSECS/SAMPLE ) 

\ INITIALISE CURRENT A/D TEMPLATE 


\ A/D INPUT TO TEMPLATES BUFFER 

\ SET TIME AXIS 


DATA. BUFFER XSECT[ 1 ] 2048 - 
SIGNAL := SIGNAL 
-.9336 ** 1324. * 

SIGNAL := SIGNAL 
. 1 SET . CUTOFF . FREQ 
SMOOTH SIGNAL := SIGNAL 
MEAN MEN := SIGNAL 
MEN - SIGNAL := 


\ CHANNEL 0 ON STACK (ANEMOMETER) 
\ DEGREES CENTIGRADE 
\ SMOOTH DATA (CYCLES, 'POINT) 

\ CENTER ON ORIGIN 


ASYST Version 3.00 

Page 1 REV. FOR 01/06/95 07:18:14.81 


QWG1NAL PAGE fS 
OF POOR QUALJ7V 


TIME SUSl 100 
XY. AUTO. PLOT 


300 


1 ] SIGNAL SUE[ 100 . 300 , 1 ] 
\ PLOT CHANNEL 0 


SIGNAL 

THAN SI := TRANS 1 ZMAG 
5 SET . # . POINTS 
1 SET . U . OPTIMA 
ZMAGO SUB! 1 , 120 , 


\ TAKE 

ZMAGO := \ MAGNI 

1 ] LOCAL. MAJIIMA 


SWAF NUM := NUM 
1 - CD / 4. / 5 
AKS [ 1 ] := AN 
DEG TRANS1 ZARG 
AKS : 3 ] := 
TRANS1 [ NUM ] 


12. / 1000. * PI * 2. 
S [ 1 3 DROP ANS [ 2 ] 
ZARGO := ZARGO [ NUM 

DEG ZARG ANS [ 3 ] : 


* 



\ 

\ 

\ 

\ 

\ 


OF FFT 


FIND INDEX AND MAX OF MAGN I T! 
INDEX OF MAX MAGNITUDE 
FREQUENCY AT MAX MAGNITUDE 

FIND ARGUMENT OF FFT 
ARGUMENT AT MAX MAGNITUDE 


DATA. BUFFER XSECTl 
SIGI := SIG1 
.04725 * 

SIGI := SIGI 


204S - 


\ SMOOTH SIGi 
MEAN MEN! : - 
TIME SUB [ 100 , 
\ XY. DATA. FLO’ 


SIGI men: 
3C0 , 1 ] 



\ channe: 


ON 


stack 




\ SMOOTH DATA ( CYCLES/PCI NT) 

ICO , 300 , 1 J 
\ FLOT CHANNEL 1 


SIGI 

FFT 

TRANS! T 
ZMAGO SUSl 
SWAF NUM : = 
1 - CD / 4. 
ANS L 4 I : 
DEG TRANS! 


RANG! ZMAG ZMAGO := 

! , 120 , 1 ] LOCAL. MAXIMA 


NUM 

/ 512. / 1000. i PI * 2 . * 
= ANS : 4 ] DROP ANS [ 5 ] 
ZARG ZARGO := ZARGO [ NUM j 


ANS [ 6 ] := 
TRANS!' [ NUM 


ZARG 


6 ] : = 


TAKE 

MAGN 



D MAX 0 
MAGN I TU 
MAX MAG 



\ FIND ARGUMENT OF FFT 
\ ARGUMENT AT MAX MAGNITUDE 


SIGI 

D I FFERENT I ATE . DATA 


.00001 + DIF! 


DIFFERENTIATE CKANNu- 


DATA . BU 

FFER : 

SIG2 := 

SIG2 

.05457 

* 

SIG2 := 

c t 

1 VJi. 


2043 


\ CHANNEL 2 ON STACK ( TKERMOCC 


SMOOTH SIG2 := 
MEAN MEN 2 := 
TIME SUBL 100 , 


SIG2 

SIG2 MEN 2 
300 , 1 ] 


\ SMOOTH DATA (CYCLES/POINT) 
- SIG2 := \ CENTER ON ORIGIN 

SIG2 SUBL 100 , 300 , 1 ] 


ASYST Version 3.00 

Page 2 REV. FOR 01/06/95 


07:13:37.93 


PAG£ f.v 

OP *00* QUAtrr 



\ XY.DATA.PLO' 


\ PLOT CHANNEL 2 

SIG2 

DIFFERENTIATE. DATA .00001 + DIF2 \ DIFFERENTIATE CHA: 


DATA. BUFFER X£ECT[ 4 j 2048 
SIG3 3IG3 
-0472E * 

SIG3 := 51G3 



SMOOTH SIG3 := SIG3 
MEAN MENS SIG3 MENS - SIG3 := 

TIME 3UBZ 100 , 300 , 1 ] SIG3 SUEZ 100 , 
\ XY. DATA. PLOT 

SIG3 

DIFFERENTIATE. DATA .00001 + DIFS 


\ SMOOTH DATA {CYCLES/POINT; 
\ CENTER ON ORIGIN 
300 , 1 ] 

\ PLOT CHANNEL 3 


\ DIFFERENTIATE CHANNEL 3 




ASYST Version 3.00 

Page 3 REV. FOR 01/06/95 07:14:37.80 


OWGINAL PAGE PS 
OF POOR QUALITY 



